from numpy import sqrt, real, imag
from matplotlib.pyplot import plot, show

def figure_2b(kappa1, alpha, beta):
	sigma_000 = sqrt(kappa1)*(alpha[:,0,0,0] + beta[:,0,0,0]);
	sigma_011 = sqrt(kappa1)*(alpha[:,0,1,1] + beta[:,0,1,1]);
	sigma_001 = sqrt(kappa1)*(alpha[:,0,0,1] + beta[:,0,0,1]);
	sigma_111 = sqrt(kappa1)*(alpha[:,1,1,1] + beta[:,1,1,1]);

	plot(real(sigma_000), imag(sigma_000),'r',linewidth=2)
	plot(real(sigma_011), imag(sigma_011),'r--',linewidth=2)
	plot(real(sigma_111), imag(sigma_111),'b',linewidth=2)
	plot(real(sigma_001), imag(sigma_001),'b--',linewidth=2); 
	show()
	pass